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Novae are thermonuclear explosions on a white dwarf surface fueled by mass 
accreted from a companion star. Current physical models posit that shocked 
expanding gas from the nova shell can produce X-ray emission but emission at 
higher energies has not been widely expected. Here, we report the Fermi Large 
Area Telescope detection of variable 7-ray (0.1-10 GeV) emission from the 
recently-detected optical nova of the symbiotic star V407 Cygni. We propose 
that the material of the nova shell interacts with the dense ambient medium 
of the red giant primary, and that particles can be accelerated effectively to 
produce 7r° decay 7-rays from proton-proton interactions. Emission involving 
inverse Compton scattering of the red giant radiation is also considered and is 
not ruled out. 

V407 Cygni (V407 Cyg) is a binary system consisting of a Mira-type pulsating red gi- 
ant (RG) with a white dwarf (WD) companion; these properties place it among the class of 
symbiotic binaries (7). Historically, although one of the more active symbiotic systems, V407 
Cyg showed an optical spectrum in quiescence dominated by the Mira-like RG (M6 III) and 
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only weak emission lines, e.g., (2). Its infrared continuum (consistent with a dusty wind) and 
maser emission (3) are detected at levels similar to other symbiotic Miras [e.g., R Aqr (4)]. 
One outstanding anomaly of V407 Cyg is a strong Li I A6707 line indicative of an overabun- 
dance of Li relative to normal Mira red giants (J, 6). Based on the 745 -day pulsation period 
of the RG (7) and the Mira period-luminosity relation (S), we adopt the distance, D = 2.7 
kpc, estimated as the mean derived from photometry in three near-infrared bands assuming an 
extinction, -Eb-v = 0.57 (2). 

A nova outburst from V407 Cyg was detected on 10 March 2010 (9); it had a magnitude 
approximately 6.9 in an unfiltered CCD image obtained at 19:08 UT. Subsequent densely sam- 
pled observations show the outburst was followed by a smooth decay, though the precise epoch 
of the nova is formally uncertain by up to 3 days due to the time gap from the pre-outburst 
image (Fig. 1). Monitoring of the source over the past two years indicates pre-outburst magni- 
tude values in the range 9-12 [see the supporting online material (SOM)]. V407 Cyg has been 
monitored optically for decades and has shown earlier signs of optical brightening on month- 
timescales by 1-2 magnitudes in the B and V bands (around 1936 and 1998) from typical V 
band magnitudes of 13-16 (2, 10,11), but the magnitude of the recent nova was unprecedented. 

Here we report on a high-energy 7-ray source (Fig. 2) positionally coincident with V407 
Cyg detected after the nova (12) during routine automated processing of all-sky monitoring data 
from the Fermi Large Area Telescope (LAT) (13). A 7-ray light curve (one-day time bins) of 
this source generated from an analysis of all LAT data reveals that the first significant detection 
(4.3a) was in fact on 10 March indicating the 7-ray activity began on the same day as the 
reported optical maximum of V407 Cyg (Fig. 1, SOM). The observed 10 March flux is up to a 
factor of 3 larger than the one-day upper limits (unless otherwise noted, 95% confidence limits 
are reported throughout) on the pre-outburst days. To further isolate the onset of detectable 
7-ray emission, we divided the 10 March data into 6-hour intervals and the first indication of 
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a signal was a marginal detection in the 12h-18h bin (2.8(t). This was followed by a highly 
significant detection (at 5.1a) in the last 6-hr bin (18h-24h) which had a peak flux that was a 
factor of ^ 3 larger than that of the marginal detection and the limits from earlier in the day 
(SOM). The initial detection of the 7-ray source by the Fermi-LAT in the latter 1/2-day of 10 
March is consistent with the time of the optical nova discovery. 

The peak flux in 7-rays (defined on one-day segments) was observed between 13-14 March, 
3-4 days after the initial 7-ray detection, and with a factor of 2 greater flux (9 x 10~^ photons 
cm~^ s~^). Analyzing data up to mid-April, the last day with a significant detection (>3cr) of 
the variable 7-ray source is 25 March, amounting to a total lifetime of activity of two weeks. 
Defining an 'active' period (14) from 10 March 18:00 to 29 March 00:00, we obtained a 7-ray 
position (Fig. 2) from the cumulative exposure that is 0.040° offset from the optical position 
of V407 Cyg, which is within the LAT 95% confidence circle (radius = 0.062°). The average 
spectral energy distribution (SED) of the Fermi-LAT 7-ray source during the defined active 
period can be described with an exponentially cut off power-law model (SOM) with a flux 
(>100 MeV) of (4.4 ± 0.4) x 10"'^ photons cm~^ s"^ (overall source significance of 18.1cr). 
A likelihood ratio test demonstrates that the addition of the exponential cut off improves the 
fit at the 4.9(7 level compared to a simpler single power-law model. We find no evidence for 
spectral variability over the duration of the active 7-ray period (SOM). Analyzing data from the 
two weeks (29 March to 12 April) following the active period collectively gives a significance 
of only 1.6(7 (flux upper limit of 0.8 x 10~^ photons cm~^ s~^), indicating that the flux has 
declined below detectability. Overall, the 7-ray source is brightest at earlier times, consistent 
with the optical behavior of V407 Cyg. The coincident localization and the observed correlated 
variability imply that the optical nova is the source of the variable 7-ray flux. 

Nuclear 7-ray lines and continuum emission from novae at ^ 1 MeV energies have been 
considered (15), but the Fermi-LAT detection of V407 Cyg shows unequivocally that novae can 
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generate high-energy (>100 MeV) 7-rays. The >100 MeV 7-ray luminosity, its spectrum, and 
light curve can be understood broadly as consequences of shock acceleration taking place in 
a nova shell. Such a nova shell is produced by thermonuclear energy release on the WD and 
initially expands freely into a very dense medium consisting of the RG wind and atmosphere 
present in the binary system prior to the nova. The radio source detected from V407 Cyg over 
22-28 March {16, 1 7), and subsequent imaging which revealed the emission to be predominantly 
extended at few milli-arcsecond resolution {18), is consistent with a picture of an extended shell 
as was found in the 2006 RS Oph nova outburst {19, 20). 

An initially spherical shell can sweep up mass from the companion RG wind asymmetrically 
and will reach a deceleration phase during which it expands adiabatically {21) with different 
temporal behavior in different directions from the WD center. Fermi acceleration of protons 
and electrons takes place in the outgoing nova shock during both the free expansion phase and 
deceleration phase, and we show that the measured 7-ray spectrum can be explained by 7r° 
decay 7-rays from proton-proton {pp) collisions or inverse Compton (IC) scattering of infrared 
photons from the RG by electrons. In both these scenarios, the 7-ray light curve in conjunction 
with the delayed rise of the X-ray flux can be explained qualitatively as a geometrical effect of 
the nova shell evolution. 

The measured optical peak magnitude w 7 of V407 Cyg (9) over one day, implies an energy 
release of > 3 x 10^^ ergs at visible frequencies (after extinction correction). The kinetic 
energy of the ejecta in the nova shell, ~ 10^'' ergs, can be estimated assuming a nova shell 
mass, Mej ~ 10~^ Mq [which is a plausible value for a massive WD, >\.25Mq {22, 25)], and 
the velocity, Vej = 3200 ± 345 km s~^, inferred from an Ha line width measurement on 14 
March (SOM). The velocity of the outgoing shock wave is initially Vsh ~ ^ej- The onset of the 
deceleration phase and subsequent evolution of the shock wave are determined by the density of 
the RG stellar wind and atmosphere surrounding the nova shell, which depends on two poorly 
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constrained parameters - the WD-RG separation (a) and the RG mass loss rate in the wind 
(Mw). As we will discuss shortly, the detection of the 7-ray flux within one day of the optical 
detection of the nova, a peak duration of 3-4 days, and subsequent decay within 15 days after 
the optical nova can be modeled as emission from the nova shell in a dense environment and 
mostly from along the WD-RG axis assuming an inverse- square law wind density profile from 
the RG center, with a ~ 10^^ cm, and ~ 3 x 10~^ Mq yr~^. Both these values are a factor 
2 smaller than those suggested previously for the V407 Cyg system (2). 

The density of particles in the RG wind is n{R) = [An {li^ + — 2aRcos 6) m]^^. 
Here, the RG wind velocity, ?s 10 km s~^, is based on optical spectra (SOM), the mean 
particle mass is m 10~^^ g, R is the distance from the WD center, and 9 is the polar angle 
relative to the WD center. The energy density in the RG radiation field is similarly, mir(-R) = 
LiK [47r(i?^ -\-a^ — 2aRcos 6) c]~^ The RG luminosity, Ljr w IO'^Lq (2, 11), is consistent with 
a spectral fit to post-nova infrared measurements with a temperature of ^ 2500 K (SOM). Near 
the WD surface {R ^ 0.01 Rq), these densities are n ~ 10^ cm~^ and mir ~ 0.01 erg cm~^, 
and increase by up to an order of magnitude when the nova shell approaches the RG surface 
(i.e., along 9 0°) at a radius trg ~ 500 Rq. An equipartition of the energy density in 
the magnetic field expected to arise from turbulent motions in the wind, to the thermal energy 
density in the RG wind with temperature w 700 K (2) gives a mean magnetic field, -Bsh(-R) = 
[327r n{R) k T^]^-^^ ~ 0.03 G in the shock wave when it is near the WD. Electrons and protons 
can be accelerated efficiently in this magnetic field (24) and interact with the surrounding RG 
wind particles and radiation. 

The time scale for pp interactions for a 7r° model to produce 7-rays in the shock wave 
is tpp Ri l/[An{R) CfTpp] ~ 2.8 x 10^ s when the nova shell is near the WD. Here, o-pp « 
3 X 10~^^ cm^ is the pp cross-section. Thus, t/tpp ~ 3% of the protons can interact to pro- 
duce 7r° emission on a time scale, t — 1 day. In an IC scenario, the cooling time scale 
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for electrons with energy Ee ^ 5 GeV that up-scatter 2500 K photons to 100 MeV is 
tic ~ (3/4) ml [aT uui{R)]-^ ~ 3.1 x 10^ s. Thus, t/tic ~ 28% of the electrons produce 
7-rays efficiently in a time scale, t — 1 day. The efficiency for 7-ray production in both the 
7r° and IC models increases substantially in the part of the nova shell that expands towards the 
RG {6 ~ 0°) and reaches the deceleration phase, by accumulating RG wind and atmospheric 
material of mass equal to Mej, at a distance ~ 7 x 10^^ cm in about 2.5 days. The efficiency 
decreases rapidly in the part of the shell that expands away from the RG {9 > 90°) because of a 
decreasing density in both the RG wind and radiation. 

In our scenario, most of the 7-rays come from the part of the nova shell approaching the 
RG. This can qualitatively explain the basic features of the 7-ray light curve (Fig. 1): (i) 
its onset within days of the optical nova and the peak flux reached in three days due to an 
increasing efficiency for pp interactions and an increasing volume of the shock-accelerated 
particles, and (ii) the decline in the flux after about five days due to weakening of the shock 
wave after reaching the deceleration phase. A highly significant one-day detection (6.5(t) of 
an increase in 7-ray flux 9 days after the nova discovery (Fig. 1) could be due to the nova 
shell hitting a part of the RG surface or a nearby remnant of high density with a size scale of 
~ 10^3 (i;^j^/i500 km s-i) (t/1 day) cm. 

A representative 7r° decay model is shown in Fig. 3 fitting the LAT data with a cosmic- 
ray proton spectrum in the form of an exponentially cut off power-law, A^p = A^p^o (Wv + 
m-p c^)~^^ e^^p/^'^p, where E^^, Wp, and rrip are the cutoff energy, kinetic energy, and mass 
of the proton, respectively. The 7-ray spectrum is well reproduced with a spectral index, Sp — 
2.15jlo;28' similar to the expected spectrum from the Fermi acceleration mechanism, with Ecp — 
32l8^ GeV (la uncertainties) (SOM). The total energy in 7-rays above 100 MeV integrated over 
the active period is a; 3.6 x 10^^ erg. The total energy in protons is Sp a; {Ep/ E^)£^, 
where ep is the mean efficiency for pp interactions and 1 / {Ep/ E^) fa 0.2 is the mean fraction 
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of proton energy transfer to 7-rays per interaction in the > 100 MeV range. The ratio of the 
total energy in protons that produce 7-rays to the kinetic energy, £p/£k ~ 9%, with ep ~ 0.2 
averaged over the 7-ray source lifetime of 15 days and the whole nova shell, similar to 1-10% 
estimated in supernova remnants (25). The ratio, however, is larger when considering dominant 
7-ray emission coming mostly from the part of the shell that expands towards the RG. 

The leptonic model is represented in Fig. 3 as the total of the IC spectrum plus a small 
contribution from bremsstrahlung emission, that arises from scattering of electrons with pro- 
tons of density n{R) in the shock wave (SOM). The exponentially cut off power-law electron 
spectrum, — A^e,o W'^" g-We/Ece^ where £Jce and are the cutoff energy and kinetic en- 
ergy of electrons (GeV), respectively, reproduces the 7-ray spectrum with Se = ~l-75lo;59 
and Ece = 3.2lo;i GeV (la uncertainties) (SOM). The total number of electrons required in 
steady state is A^e,o ~ 4 x 10^^, with a mean energy of 8.7 GeV which is larger than Ece be- 
cause of the steep spectrum. The total energy in electrons over 15 days of 7-ray emission is 

~ 4 X 10^^ ergs, averaged over the nova shell. Thus the total energy in electrons that produce 
7-rays in the IC model is a small fraction of the kinetic energy in the shell {E^/S-k ~ 0.4%). 

X-ray emission detected from V407 Cyg with the Swift X-ray Telescope (XRT) as early as 
three days after the onset of the optical nova (Fig. 1) is likely due to shock-heating of ambient 
gas {26,27). The X-ray flux starts rising significantly about two weeks after the nova, coinciding 
with when the 7-ray flux declines below detectability. In our geometric scenario, the sharply 
rising X-ray flux is due to the increasing volume of shocked gas in the nova shell expanding in 
the direction away from the RG. The X-ray flux peaks about 30 days after the explosion, and its 
subsequent slow decline is consistent with the longer time scale of the deceleration phase. 

The F(?rm/-LAT detection of V407 Cyg was a surprise, and adds novae as a source class 
to the high-energy 7-ray sky. The particle acceleration mechanism and the 7-ray emission 
scenarios outlined here require the mass donor to be a red giant, i.e., a nova in a symbiotic 
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system. Interestingly, several symbiotic stars are known to be recurrent novae (i.e., systems 
observed to have undergone multiple thermonuclear runaways within the last century), and 
recurrent novae are often considered candidate progenitors of Type la supemovae (28). V407 
Cyg may also belong to this class of binaries, and we have adopted parameters that are consistent 
with such a classification in modeling the 7-ray emission. These sources can in general have 
dramatic influence on the local interstellar medium and Galactic cosmic rays but few binary 
systems with a WD are known to have a similar environment, hence we expect 7-ray novae to 
be rare. 

References and Notes 

1. S. J. Kenyon, in Evolutionary Processes in Interacting Binary Stars, Edited by Y. Kondo et 
al, lAU Symp. 151, 137 (1992). 

2. U. Munari, R. Margoni, R. Stagni, Mon. Not. R. Astron. Soc. 242, 653 (1990). 

3. S. Deguchi, J. -I. Nakashima, T. Miyata, Y. Ita, Publ. Astron. Soc. Japan 57, 933 (2005). 

4. R. J. Ivison, E. R. Seaquist, R J. Hall, Mon. Not. R. Astron. Soc. 269, 218 (1994). 

5. S. Y. Shugarov, A. A. Tatamikova, E. A. Kolotilov, V. I. Shenavrin, B. F. Yudin, Baltic 
Astronomy, 16, 23 (2007). 

6. A. A. Tatarnikova et al., Mon. Not. R. Astron. Soc. 344 1233 (2003). 

7. L. Meinunger, Mitt. Veranderl. Sterne. 87, 111 (1966). 

8. I. S. Glass, M. W. Feast, Mon. Not. R. Astron. Soc. 199, 245 (1982). 

9. K. Nishiyama, F. Kabashima, lAU Central Bureau for Astronomical Telegrams, reported 
by H. Maehara, No. 2199 (2010). 

8 



10. C. Hoffmeister, Veroff. Stemw. Sonneberg. 1, 295 (1949). 

11. E. A. Kolotilov, V. I. Shenavrin, S. Y. Shugarov, B. F. Yudin, Astron. Rep. 47, 777 (2003). 

12. C. C. Cheung et al. (Fermi-LAT collaboration). The Astronomer's Telegram 2487, 1 (2010). 

13. W. B. Atwood etal. (Fermi-hAT collahovation), Astrophys. J. 697, 1071 (2009). 

14. The window was extended by three days to include any late low-level emission which would 
not be detectable on any individual daily bin. 

15. M. Hemanz, in Classical Novae, 2nd ed.. Edited by M. F. Bode & A. Evans, Cambridge 
University Press (2008), p. 252. 

16. I. Nestoras etal.. The Astronomer's Telegram 2506, 1 (2010). 

17. G. C. Bower, R. Forster, C. C. Cheung, The Astronomer's Telegram 2529, 1 (2010). 

18. M. Giroletti etal. The Astronomer's Telegram 2536, 1 (2010). 

19. T. J. O'Brien et al.. Nature 442, 279 (2006). 

20. M. P. Rupen, A. J. Mioduszewski, J. L. Sokoloski, Astrophys. J. 688, 559 (2008). 

21. L. 1. Sedov, Similarity and Dimensional Methods in Mechanics, New York: Academic 
Press (1959). 

22. S. Starrfield, W. M. Sparks, J. W. Truran, M. C. Wiescher, Astrophys. J. Supp. 127, 485 
(2000). 

23. O. Yaron, D. Prialnik, M. M. Shara, A. Ysovetz, Astrophys. J. 623, 398 (2005). 

24. V. Tatischeff, M. Hemanz, Astrophys. J. Lett. 663, LlOl (2007). 

9 



25. A. A. Abdo et al. (Fermi-LAT collaboration), Astrophys. J. Lett. 710, L92 (2010). 

26. J. L. Sokoloski, G. J. M. Luna, K. Mukai, S. J. Kenyon, Nature 442, 276 (2006). 

27. J. J. Drake et al. Astrophys. J. 691, 418 (2009). 

28. W. M. Wood-Vasey, J. L. S,o\io\o^\±, Astrophys. J. Lett. 645, 53 (2006). 

29. A. A. Abdo et al. (Fermi-LAT collaboration), Astrophys. J. Supp. 188, 405 (2010). 

30. The Fermi LAT Collaboration acknowledges support from a number of agencies and insti- 
tutes for both development and the operation of the LAT as well as scientific data analysis. 
These include NASA and DOE in the United States, CEA/Irfu and IN2P3/CNRS in France, 
ASI and INFN in Italy, MEXT, KEK, and JAXA in Japan, and the K. A. Wallenberg Foun- 
dation, the Swedish Research Council and the National Space Board in Sweden. Additional 
support from INAF in Italy and CNES in France for science analysis during the operations 
phase is also gratefully acknowledged. P. K. acknowledges support from Al ASCR (On- 
drejov) in the Czech Republic. We acknowledge with thanks the variable star observations 
from the AAVSO International Database contributed by observers worldwide and used in 
this research. 



10 



_ 12 
6- 

2 4f 



£ 8: 
1 9- 

10 ; 

11 : 

b 20 

15 
5^ lOl: 

DC 5 





LAT 

>100 MeV 



-*I-Lh 1 1 1 — I 1 1 1 1 — I 1 1 1 1 — I 1 1 1 H 



V - 



H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 h 



XRT 
0.3-10 keV 



♦ • 



10 20 30 

Days since 10 March 2010 

Fig. 1. Light curves of V407 Cyg in 7-rays from the Fermi-hAI (top), optical (middle), and 
X-rays from Swift (bottom). Vertical bars indicate la statistics errors for all data (the errors are 
smaller than the points in the optical). For the 7-ray data, gray arrows indicate 2a upper limits, 
and horizontal bars indicate the one-day binning. In the optical, unfiltered (black), Rq (red), 
and V (green) band magnitudes are shown (SOM). The vertical dashed blue line indicates the 
epoch of the optical nova detection; the 7-ray peak occurred 3-4 days later. 
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Fig. 2. Fermi-LKT 7-ray counts (200 MeV-100 GeV) map from 10 March 18:00-29 March 
00:00 2010 in the region around the position of the 7-ray nova source V407 Cyg (marked by 
black cross) at / = 86.958°, b = -0.513° (R. A. = 315.551°, Dec. = +45.737°, J2000.0). The 
map was adaptively smoothed by imposing a minimum signal-to-noise ratio of 7. The closest 
known 7-ray source is contained in the 1st year LAT catalog (IFGL J21 11.3-1-4607; marked 
by the black circle) (29), ~1.5° away from the star's optical position. The bright source at 
(/ = 78.2°, b = +2.1°) is LAT PSR J2021-I-4026. 
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Fig. 3. SED of V407 Cyg in MeV/GeV 7-rays measured by the Fermz-LAT over the period 10 
March 18:00 - 29 March 00:00 2010. Vertical bars indicate la statistical errors, arrows indicate 
2a upper limits, and horizontal bars indicate energy ranges. The best-fit vr" (black solid line) 
and leptonic (blue dashed line) models are indicated. 
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Supporting Online Material (SOM) 

Materials and Methods: Fermi-L AT Data Analysis 

The analysis of the LAT data was performed using the Fermi Science Tools v9rl5 package 
available from the Fermi Science Support Center (FSSCjj. The standard onboard filtering, event 
reconstruction, and classification were applied to the data {SI), and the high-quality ("Pass 6 dif- 
fuse") event class is used. Throughout the analysis, the "Pass 6 v3 Diffuse" (P6_V3_DIFFUSE) 
instrument response functions (IRFs) are applied. 

Events in the range 0.1-100 GeV were extracted from a 10° x 10° square region of interest 
(ROI) centered on the known location of V407 Cyg. To greatly reduce contamination from the 
Earth albedo photons, we excluded time periods when the 10° region around V407 Cyg was 
observed at a zenith angle greater than 105° and for observatory rocking angles of greater than 
52° for observations after mission elapsed time (MET) 273628805 (rocking angles of greater 
than 43° are excluded for observations prior to this time). 

The gt like likelihood fitting tool was used throughout to perform a binned spectral anal- 
ysis, wherein a spectral-spatial model containing point and diffuse sources is created and the 
parameters obtained from a simultaneous maximum likelihood fit to the data. The model was 
constructed by including the 5 brightest point sources from the IFGL catalog (52) within 15° of 
the center of the ROI: IFGL J202 1.0-1-3651, IFGL J202 1.5-1-4026, IFGL J2032.2-I-4127, IFGL 
J2030.0-I-3641, and IFGL J21 11.3-1-4607. The other IFGL sources within the ROI were not 
considered as they are intrinsically faint and not detectable on the timescales of this outburst. 
The standard models for the Galactic diffuse emission {gll jiemjj 02 . f it) and isotropic back- 
groundo currently recommended by the LAT team were also incorporated into the model. The 
first three IFGL sources listed above are known pulsars and were modeled by an exponen- 



'See the FSSC website for details of the Science Tools: http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
^Descriptions of the models are available from the FSSC: http://fermi.gsfc.nasa.gov/, 



tially cut off power-law, while the remaining two sources were modeled by a single power-law 
with parameters initially set to the values obtained in the IFGL catalog (52). In the fitting, the 
Galactic diffuse emission model was scaled by a single power-law with free normalization and 
index in order to allow for small spectral errors in the model of diffuse emission. All source 
parameters were left free when investigating the average spectral behavior over the duration of 
the entire outburst. The parameters of the background point sources were then fixed to their 
average values when performing likelihood analysis on shorter timescales. 

The average spectral behavior 

A detailed re-analysis of the LAT data on the reported peak flare days (13 and 14 March) (S3) 
confirms a new 7-ray source not previously reported by the LAT. Fitting the source with a single 
power-law spectrum yields a flux (>100 MeV) = (1.3±0.2) x 10^^ photons cm~^ s~^, and slope, 
r = 2.2 ± 0.1, with a test statistic (TS, (S4)) of 111, where the source significance, ~ y/TS = 
10.5(j. Subsequent analysis indicated that the source was still detectable on timescales of a day 
through 25 March. Consequently, to investigate the average spectrum and maximize the statis- 
tical significance of the detection, we define a time window extended back to the approximate 
onset of the optical outburst seen in V407 Cyg (10 March 18:00). To include any late low-level 
emission which would not be detectable on an individual daily basis, the window was extended 
out to 29 March 00:00. This gives a MET range of 289936803-291513602 and the results of 
using this defined "active period" are described below. 

Using the gtf indsrc tool on this active period gives an improved localization over the 
initially reported one (72) at (J2000.0) R.A. = 315.551°, Dec. = 45.737° (I = 86.958°, b = 
-0.513°) with a 95% error radius of 0.062°. This 7-ray position is 0.040° offset from the 
optical position of V407 Cyg; throughout the remainder of the analysis the nominal optical 
position of V407 Cyg is used. The previously discussed model was used with the gt like tool 
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to identify the average spectral properties of the source during this period. The source was best 

fit by an exponentially cut off power-law model (dN/dE oc E^^e^'^^^^"'^). The flux (>100 
MeV) obtained is (4.4 ± 0.4 (stat) ± 0.2 (syst)) x 10"'^ photons cm"^ s"\ with a photon index, 
r = 1.5 ± 0.2 (stat) ± 0.04 (syst), and a cutoff energy, E^ = 2.2 ± 0.8 (stat) ± 0.2 (syst) GeV; 
see below for discussion of how the systematic errors are estimated. The source is detected 
with a TS = 326.9 (~18.1cr). The cutoff power-law was compared to a single power-law model 
by a likelihood ratio test. This gives TS = — 2Alog(Likelihood) = 23.6 and indicates that the 
addition of the exponential cut off improves the fit at the 4.9a level compared to a single power- 
law model. 

A number of effects are expected to contribute to the systematic errors. Primarily, these 
are uncertainties in the effective area and energy response of the LAT as well as background 
contamination. These are currently estimated by using outlier IRFs that bracket the nominal 
ones in effective area. These are defined by envelopes above and below the P6_V3 -DIFFUSE 
IRFs by linearly connecting differences of (10%, 5%, 2%) at log(E/MeV) of (2, 2.75, 4), 
respectively. 

Calculating the upper limit 

Using the ~19 months of all-sky monitoring data prior to the onset of activity from this source 
described above it is possible to calculate an average flux upper limit for V407 Cyg in this 
time range. To this end a new source model is constructed which comprises all 38 IFGL point 
sources (52) within 15° of the center of the ROI and the standard models for the isotropic 
and Galactic diffuse emissions. All of the IFGL sources were modeled with single power-law 
spectra except for those known to be pulsars in which case an exponentially cut off power- 
law was applied. An additional point source is inserted at the location of V407 Cyg with an 
exponentially cut off power-law fixed to the average spectral parameters achieved in the fit to its 
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active period. The normalization of the source is allowed to be free and a maximum likelihood 
fit performed. Applying the method of (55) the 95% upper limit on the flux (>100 MeV) is 
calculated to be 6x 10~^ photons cm~^ s~^. 

Exploring the evolution of the outburst 

To investigate the evolution of the source flux the data were extracted in one-day segments. A 
likelihood analysis was performed on each of these segments using the optimized background 
source model and an exponentially cut off power-law spectrum for the source of interest. The 
spectral model parameters were fixed to the best obtained average values and only the source 
normalization and the diffuse background sources were allowed to vary. In each segment the 
source TS, flux (> 100 MeV), and 95% flux (> 100 MeV) upper limit were calculated (Table SI). 
This analysis indicated that the first day in which a significant detection for the source was 
achieved on March 10 and that the source was detected (TS>9) on daily timescales up to and 
including 25 March. To further explore the onset of detectable 7-ray emission around the date 
of the optical detection of the nova (10 March), we divided the data into 6-hr bins. This analysis 
indicated that the 7-ray signal was isolated to the latter half of this day (Table S2). 

Search for spectral variability 

We split the Fermi-LAI observation period into two time segments with roughly equivalent 
statistical significance (TS value) for the detection of the source: 10 March 18:00 - 14 March 
12:00 (a) and 14 March 12:00 - 29 March 00:00 (b). The best fit spectral index and cutoff 
energy of the power-law with exponential cutoff are F = 1.6 ± 0.2 and Ec = (3.0 ± 0.9) 
GeV for the first segment (a), and T = 1.3 ± 0.3 and = (1.7 ± 0.6) GeV for the second 
segment (b) (Fig. SI). The flux (>100 MeV) varies from (7.2 ± 3.1) x 10"'^ photons s"^ cm-^ 
to (3.3±0.8) xlO"'' photons s~^ cm~^ between the two periods while the index and cutoff 



4 



energy are statistically equal. For comparison, we fitted the total flux of the power law with 
exponential cutoff with the segment (b) data by keeping the index and cutoff energy fixed to the 
values obtained with the best fit of segment (a). It resulted in a difference in log(likelihood) of 
0.93. This difference corresponds to a significance of lAa which means that we did not detect 
any spectral variability. 

Parameters of the 7-ray Emission Models 

We calculate the 7r° production by pp interactions following the prescription of (S6). We use 
cosmic-ray proton spectra of the form, A^p = A^p o {Wp + rUp c^)~*p e^^p/^-^p (proton GeV~^), 
where Wp is the kinetic energy of protons (GeV) and rrip the proton mass. We fit the normal- 
ization (A^p,o), spectral index (sp), and the cutoff energy (Ecp), with the Fermi-LNI data. The 
7r° emissivity was calculated assuming a solar metallicity. We apply a corresponding nuclear 
enhancement factor of en = 1-84 to the 7r° emissivity (57). 

The resulting best fit spectral model presented in the main text has a spectral index, Sp — 
2.15lo;28 ^i^d ^cp — 32^8^ GeV (lu uncertainties). Fig. S2 shows the confidence region of 
the spectral index and cutoff energy fit to the Fermi-LAT data. Although the uncertainty of 
the spectral index is large, its best fit value corresponds to the canonical slope of cosmic-ray 
protons (sp ~ 2) in the first-order Fermi acceleration process. Note that a cutoff energy larger 
than ~100 GeV is not excluded for spectral indices larger than ~2.5 (95% confidence level). 

The observed 7-ray flux in the 7r° model can be reproduced with a total number of cosmic- 
ray protons, JNpdWp = 2.0tJ:^ x 10^^ (^h/4 x 10^ cm-^^)-^ (D/2.7kpc)2 in a steady state. 
The total energy in protons is / Wp Np dWp = 6.9ty x 10^^ (^^/4 x IQ^ cm-^)"! (D/2.7 kpc)^ 
ergs. Here, nn ~ 4 n{R), is the density of target particles in hydrogen gas in the shock region 
and D^2.7 kpc is the adopted distance to V407 Cyg (S8). 

The inverse Compton and bremsstrahlung emissions in the leptonic model are calculated 
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using the method presented in (S9). In our steady state calculation, cosmic-ray electrons in the 
nova shell interact with infrared photons from the red giant (RG; modeled as a black-body with 
a temperature of 2500 K and a radius of 500 Rq) at a distance of ~ 10^^ cm from the RG. At this 
distance, the target ion density for the bremsstrahlung process is ~ 4 x 10^ cm~^ as in the 
7r° model. We use cosmic-ray electron spectra of the form, A^'c = Nefl e^^'=/^='' (electron 
GeV~^), where We is the kinetic energy of electrons (GeV). The normalization factor (A^e,o)> 
spectral index (se), and cutoff energy (Ece) are varied to fit the Fermi-LAT data. 

The resulting best fit leptonic spectral model presented in the main text (Fig. 3) has a spec- 
tral index, = -l.TbtlH and = 3.21^;? GeV (la uncertainties) (Fig. S3). Fig. S4 
shows the confidence region of the spectral index and cutoff energy fit to the Fermi-LAT data. 
It presents two optimal zones (confidence level 68%), around Se = —1-75 and Se = 0.25. 
The best fit electron level leads to a total number of cosmic-ray electrons in a steady state of 
JNedWe = 0.391:I%I X 10^=^ (D/2.7kpc)2 and the total energy in electrons of J W^N^dWe 
= 0.54lo1i3 X 10"^^ (D/2.7kpc)^ ergs. Note that the uncertainty in the electron spectral index 
is rather large, and the canonical slope (se ~ 2) expected to arise from the first-order Fermi 
acceleration process cannot be excluded (99% confidence level). 

The difference in the log(likelihood) value for the best fit leptonic model and that of the best 
fit 7r° model is ~0.05. This difference corresponds to a significance of ~0.3cr and neither model 
is statistically preferred over the other. 

Optical and Infrared Photometry of V407 Cyg 

The optical nova of V407 Cyg was discovered by Nishiyama and Kabashima (SIO) using an 
unfiltered CCD image obtained on 10.797 March 2010 UT at the Miyaki Argenteus Observatory 
in Japan with a 105-mm f/4.0 lens telescope. The pre-outburst image from 7.859 March shows 
the source 2.5 magnitudes fainter. There are uncertainties in the actual epoch of the nova due to 
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the three-day gap in the observations. Previous images obtained with the same equipment dating 
back to April 2008 (Fig. S5) show slow light variations with the total amplitude ~3 mag, which 
may be attributed to a combination of the Mira-type pulsations of the red giant and activity of 
the white dwarf. 

The bright outburst triggered follow-up photometric observations with a range of instru- 
ments and the observations shown in the main paper are summarized here (see Table S3). 
Many CCD measurements were contributed by enthusiasts from the American Association of 
the Variable Star Observers (AAVSO), and those who contributed V and Rc band CCD pho- 
tometry presented here were: G. Belcheva (Bulgaria), S. Dvorak, M. Halderman, G. Sjoberg, 
D. Trowbridge (USA), T. Kantola, M. Luostarinen, A. Oksanen, J. Virtanen (Finland), D. Lane 
(Canada), S. O'Connor (Bermuda), S. Padovan (Spain), and A. Smirnov (Russia). Observa- 
tions from the Kwasan Observatory and KANATA telescope are also reported. Small system- 
atic discrepancies between CCD measurements obtained with different telescopes (resulting 
from different choices of comparison stars and mismatches between instrumental and standard 
bandpass) were compensated when necessary by using the well sampled AAVSO Bright Star 
Monitor (BSM) light curve as the primary data set. The BSM photometry used an ensemble of 
comparison stars in the 8-9 mag range calibrated using the Landolt standards (Sll). 

Fig. S5 shows the combined optical light curve in V and Rc bands and the unfiltered CCD 
measurements by K. Nishiyama and F. Kabashima. A subset of these data are shown in the 
main paper. The magnitude zero point for the unfiltered measurements was calculated assuming 
Rc = 6.571 for the comparison star HIP 103871 as indicated by the BSM photometry, resulting 
in the peak magnitude ~0.5 mag brighter than the value reported in the discovery (SIO). The V- 
band light curve peak of 7.8 mag on 11 March 2010 20:28 communicates a larger brightening, 
in comparison with the (sparsely covered) historical magnitudes of 14-15, than is apparent in 
the unfiltered photometric data. 
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V407 Cyg was also observed in the infrared nine and twenty-five days after the nova ex- 
plosion using the InSb-photometer (512) at the 1.25-m telescope of the Crimean Laboratory 
(Sternberg Astronomical Institute). The JHKLM band magnitudes observed are typical for 
this system around maximum of the Mira pulsation cycle (although with slightly bluer colors), 
indicating that the RG was still dominating in the infrared (SI 3). Assuming the Galactic extinc- 
tion (-Eb-v = 0.57) adopted in this paper (58), we converted magnitudes to flux densities (S14) 
and fitted the infrared spectra with Planck's law, which gave a RG temperature of 2500 K at 
both epochs. 

Optical Spectroscopic Measurements 

The ejecta velocity, Vej = 3200 ± 345 km s~^, quoted in the main paper was taken as the half 
width zero intensity of the Ha line measured in a spectrum obtained at the Castanet-Tolosan 
Observatory and made available online by C. Buij^. The measurement was done by fitting a 6th 

o 

order polynomial to the continuum only in the range A6420-6655A, and subtracting from the 
observed spectrum. 

Optical spectra of V407 Cyg obtained with the Nordic Optical Telescope (program P40- 
423) with resolution 2.2 km s^ (resolving power = 67000) showed that the O I] A6363A line 
displayed a narrow emission peak at —54 km s^^ (Fig. S6) on 31 March with FWHM = 26 km 
s^ , as well as strong absorption on Na I AA5889, 5895 A (components at —58 and —51 km 

o 

s^^). Ca I A4226.7A shows a single emission component, centered at —55 kms^^, with FWHM 
of 7 km s^^ that is clearly from the stellar wind. The Balmer lines show strong, broad P 
Cyg profiles centered at —60 km s~^, indicating absorption through the wind of the RG. The 
permitted and forbidden Fe-peak lines (i.e., Fe II) display narrow components and the strongest 
permitted lines display P Cyg profiles with absorption component displacements indicating a 
■'http;//astrosurf.com/aras/V407Cyg/v407cyg.htm 
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RG wind velocity, I'w ~ 10 km s ^. 
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Fig. SI. Spectral energy distribution of V407 Cyg in MeV/GeV 7-rays measured by the Fermi- 
LAT over the period from 10 March 18:00 to 14 March 12:00 (a) and 14 March 12:00 to 29 
March 00:00 2010 (b). Horizontal bars indicate energy ranges, vertical bars indicate la statis- 
tical errors, and arrows indicate 2a upper limits. The best fit of a phenomenological power-law 
with exponential cutoff (dashed line) is presented for the two periods. 
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Fig. S2. Confidence region of the spectral slope (Sp) and logarithm of the cutoff energy (Ecp in 
GeV) fit of the Fermi-LAT data for V407 Cyg with the cosmic-ray proton spectrum. The star 
indicates the best-fit values. 
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Fig. S3. SED of V407 Cyg in MeV/GeV 7-rays measured by the Fermz-LAT over the period 10 
March 18:00 - 29 March 00:00 2010. Vertical bars indicate la statistical errors, arrows indicate 
2(7 upper limits, and horizontal bars indicate energy ranges. The best-fit leptonic model with the 
separate contributions from the IC (dot-dashed line) and bremsstrahlung (dotted line) spectra 
indicated, as well as their total (dashed line). 
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Fig. S4. Confidence region of the spectral slope (Sg) and logarithm of the cutoff energy (Ece in 
GeV) fit of the Fermi-LAT data for V407 Cyg with the cosmic-ray electron spectrum. The star 
indicates the best-fit values. 
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Fig. S5. Optical light curve of V407 Cyg extending back to April 2008. The blue vertical line 
indicates the epoch of the nova discovery on 10 March 2010. 
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Fig. S6. The O I] A6303A line in the heliocentric radial velocity observed in V407 Cyg with 
the Nordic Optical Telescope on 31 March 2010. 
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Days since 


TS 


Flux (>100 MeV) 


10 March 2010 




[lO"'^ photons cni~^ s~^] 


^.5 


1.5 


<4.4 


-3.5 


0.0 


<1 5 


-2.5 


0.0 


<2.2 


-1.5 


2.2 


<3.6 


-0.5 


0.0 


<1.4 


0.5 


18.9 


4.3 ± 1.6 


1.5 


56.9 


7.2 ± 1.7 


2.5 


12.0 


3.7 ± 1.5 


3.5 


68.4 


9.0 ± 1.9 


4.5 


47.8 


8.7 ± 2.0 


5.5 


27.8 


7.3 ± 2.0 


6.5 


16.7 


3.3 ± 1.3 


7.5 


9.7 


3.7 ± 1.8 


8.5 


11.7 


2.7 ± 1.2 


9.5 


42.2 


7.6 ± 2.0 


10.5 


11.6 


2.6 ± 1.2 


11.5 


15.9 


3.7 ± 1.6 


12.5 


0.5 


<2.8 


13.5 


7.4 


<6 9 


14.5 


19.7 


4.9 ± 1.7 


15.5 


13.3 


4.9 ± 1.8 


16.5 


0.0 


<3.5 


17.5 


4.8 


^5 7 


18.5 


0.0 


<\ 8 


19.5 


2.4 


<3.5 


20.5 


0.2 


^3 2 


21.5 


0.9 


<3.7 


22.5 


0.2 


<4.2 


23.5 


0.0 


^3 1 

\ J . J- 


24.5 


0.0 


<2.0 


25.5 


2.1 


<A 3 


26.5 


0.0 


<3.3 


27.5 


1.9 


<3.2 


28.5 


0.0 


<\ 8 


29.5 


0.0 


<2.8 


30.5 


0.2 


<3.2 


31.5 


0.0 


<1.7 


32.5 


0.3 


<2.2 


33.5 


2.1 


<4.7 


34.5 


0.0 


<2.2 


35.5 


3.1 


<5.5 


36.5 


0.0 


<1.9 


37.5 


1.6 


<3.5 


38.5 


0.2 


<2.9 



Table SI. Daily LAT 7-ray test statistic (TS) and flux values (95% confidence upper limits when TS<9) presented 
in the main paper. The dates indicated are the centers of the one-day bins. 
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Time 


TS 


Flux (> 100 MeV) 


Interval 




[10"'^ photons cm~^ s~^] 


00h-06h 


0.0 


<3.9 


06h-12h 


0.0 


<6.4 


12h-18h 


8.0 


6.2 ± 3.7 


18h-24h 


32.0 


17.4 ±5.3 



Table S2. Six-hour 7-ray test statistic (TS) and flux values (or 95% confidence upper limits) for 
10 March. Note that a flux is reported for the 12h-18h bin despite showing a TS<9 as adopted 
throughout; the corresponding upper limit was <15.1 x 10~^ photons cm~^ s~^. 



Telescope 


Filter(s) 


Camera 


Astrokolkhoz Observatory AAVSO BSM 60-mm, USA 


VRc 


SBIG ST-8XME 


Kwasan Observatory 250-mm, Japan 


VRc 


SBIG ST-7XME 


Hiroshima University KANATA 1.5-m, Japan 


V 


TRISPEC (S15) 


Miyaki Argenteus Observatory l()5-mm, Japan 


unfiltered 


SBIG STL6303E 



Table S3. Instruments used for the optical photometric monitoring observations (Fig. S5). 
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